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Abstract. We introduce a generalization of Blum-Shub-Smale machines on the standard real numbers 
R that is allowed to run for a transfinite ordinal number of steps before terminating. At limit times, 
register contents are set to the ordinary limit of previous register contents in R. It is shown that 
each such machine halts before time w” or diverges. We undertake first steps towards estimating the 
computational strength of these new machines. 


1 Introduction 


In the spirit of ordinal computability - the study of classical models of computations generalized to 
transfinite ordinal numbers - we study a variation of the Blum-Shub-Smale (BSS) machine introduced 
in [1]. In contrast to established models of ordinal computability, such as Hamkins’ and Lewis’ infinite 
time Turing machines (ITTMs) [3] and Koepke’s ordinal Turing machines (OTMs) [4], these machines 
employ real numbers in the classical continuum R as opposed to elements of Baire space “w or Cantor 
space “2. The topological differences matter as soon as we consider limits (see below). Variations thereof, 
be it in allowing infinitely many registers or changes in the limit behavior, might very well change the 
computational strength. In this paper, we aim for the “weakest” possible generalization of BSS machines 
into ordinal time. We believe that already this restricted model shows interesting properties. One could 
consider many variations and strengthenings, which will also affect the computational strength. 


Our machines have a finite number n of registers, each containing a real number. Generalizations to 
other fields and rings are possible but shall not be of concern to this paper. The computation is steered 


by a finite program P Cw x {f | f : R” a R”} x {0,1} x w x w, containing commands of the 
form (i, ¢,j,k,1), where i is the index of the command at hand, ¢ is a rational map (with rational 
coefficients), and j tells us if the command represents a computation node or a branch node. In case 
j = 0, we are at a computation node, the register content x € R” is replaced by ¢(x) and the next 
command (index i + 1) is carried out next. The values of k and l are ignored in this case. Otherwise, 
j = 1 and we are at a branching node. This means that the register content is left unchanged and, 
depending on whether (x) > 0, the next command will be the one with index k. If on the other hand 
(x) < 0, command number l is carried out next. We can assume the indices of a given program’s 
commands to form an initial segment of the natural numbers and that no index appears twice. In case 
a command index is called for which no command in the program exists, the computation halts. 


Note 1. As a minor technical detail we would like to note that, as in the original paper [1], we avoid 
discontinuity points of rational functions by putting decision nodes before each computation or decision 
node to check whether the denominator of the rational function to be evaluated is 0. If not, we continue 
as planned, if yes, an infinite loop is entered and the computation diverges. 


So far, we have outlined a standard BSS machine with the additional restriction that the rational 
functions present in computation and branching nodes do not allow for arbitrary real coefficients. 
We add irrational coefficients in form of parameters later on. We now make our machines access the 
transfinite: In order for the machine to run for infinitely many steps, we have to define the register 
content at limit times. In the established theories of infinite time or ordinal Turing and register machines, 
often an inferior or superior limit is used for this purpose. Instead, we want to restrict ourselves here 


to ordinary limits of sequences of real numbers. This immediately implies that there will be situations 
where an infinite time BSS machine will, e.g., be properly defined at any finite time but not at the first 
limit time w because the register contents do not converge. We can imagine the machine to “crash” in 
such a case and say that for such a combination of program and input no valid computations exists. 
In case of converging register contents we also have to come up with a command that is carried out at 
a limit time. For this we will indeed use the inferior limit, i.e. the command with the least index that 
was used cofinally often below the limit. Note that we do not introduce a dedicated limit state. 


Let us put things together in the following definition. 


Definition 1. Letn € w be a number of registers. Let k < n and let P be an+k-register BSS program. 

The infinite time BSS machine computation (ITBM computation) with parameters pi, p2,...,px% E€ R 

by P on some input x € R” is defined as the transfinite sequence (Ct)reo = (R(t), I(t))reo € °(R"** xw) 

where: 

(a) 0 € Ord or 6 = Ord 

(b) R(O) = (a, pi, p2,..., pk) and I(0) = 0 

(c) (computation node) Ift < 6 and I(t) = i let (i, ġ,0, k,l) be the command of P with index i. Then 
R(t+1) = ¢(R(t)) and (+1) =i4+1. 

(d) (branching node) Ift < 0 and I(t) = i let (i,¢,1,k,1) be the command of P with index i. Then 
R(t+1) = R(t) and if furthermore ¢(R(t)) > 0 then I(t +1) = k; if on the other hand $(R(t)) < 0, 
then I(t +1) =1. 

(e) Ift <@ is a limit and y = lims4+R(s), then R(t) = y and I(t) = lim infs_; I(s). 

(f) If 0 < Ord, then @ is a successor ordinal and I(0 — 1) calls a command index that is not in P (the 
machine halts). 


We define ITBM computable functions on the reals: 


Definition 2. A function f : domf > Y where domf,Y C R is called ITBM computable in param- 
eters p1, p2,..., Pk if there is an at least k + 1-register ITBM program P s.t. for every x € domf the 
computation by P on input (x,0,0,...,0,p1,p2,.--,; pk) halts and the final register content is of the form 
(f(x),-,°;---57). On input x ¢ domf the computation is required to diverge. We call such a function 
ITBM computable if k = 0, i.e. no real parameters are necessary. 


The use of one limit step enables us to compute the classical elementary functions that are defined by 
power series as illustrated by the following examples. While such functions as the exponential function 
can be computed in classical recursive analysis, they are not computable in the standard BSS model 


[2]. 


Example 1. The exponential function e : R > R, x > e” = ko = is ITBM computable: Define a 5- 
register program that computes the desired function if |x| < 1. Later we shall describe the modifications 
necessary to work for any zx. 


Algorithm 1 


input Ry = z; 

set Ro := Ra := Ra := Rs := 1; (initialize) 
call loop; 

loop: 


if Ro = 0 then set Ri := Rs and halt, else continue; 
set Ra := gy; (store 1, where i is the current iteration) 
set Rg := Rs x Ra; (store +) 
set Ro := Rə x Ri; (store x’) 
os 
set Rs := R (store 4) 
call loop; 


If |x| < 1, all register contents converge and Rs contains the desired output at time w, which is correctly 


recognized when R2 = 0. We can adapt the algorithm for |x| > 1 by adding a case distinction in the 


beginning and, in case |x| > 1, save = in register three. Then register three converges also at limit 
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times. Of course, the command that updates Rs inside the loop has to be changed accordingly. 


gk 


GEJ! and the cosine function cos : R > 


Example 2. The sine function sin: R > R, x + J k-o( 1)* 


R, s > oP o( 1)” y are ITBM computable. This is proven by the previous example and the fact 


that (—1)* can be recovered from (—4)* and (4)*, both of the latter which can be convergently stored 
and updated in separate registers. 


Note 2. Common to this type of examples is that some tricks are necessary to make all registers 
converge at limits. Auxiliary registers used for scratch work often do not contain converging content. If 
their content is bounded, however, one can simply divide the register content cofinally often by a fixed 
number and keep track of how often this division has occurred. Unbounded contents are best stored as 
their multiplicative inverse. Both approaches can be imagined as pushing the relevant data contained in 
a register into increasingly later places in their decimal/binary expansion. Compound limits like w-w are 
an additional problem, as scratch registers cannot be set to arbitrary values after limit times without 
sacrificing convergence at the compound limit. However, in every limit stage towards a compound limit 
the register content will be bounded if treated like above. So, in order to ensure convergence at the 
compound limit, these bounds themselves need to converge. See the next chapter for an example. 


2 Clockable ordinals 


Since Hamkins’ and Lewis’ paper on ITTMs [3], determining those ordinals that appear as halting times 
on empty inputs has proved to be important for the study of machine models. Since our machines do 
not halt at limit times, we are interested in machines that run for some limit number a many steps 
and halt in the next step: 


Definition 3. An ordinal a is called ITBM clockable if there is an n E€ w and an n-register ITBM 
program that halts on input 0 E€ R” in exactly a+1 many steps. 


The algorithms above prove that w is clockable, but let us establish this anew with an algorithm that 
uses only one register. 


Lemma 1. The first limit ordinal w is ITBM clockable. 
Proof (by algorithm). 


Algorithm 2 


set Ry := 1; 
call loop; 


loop : 
if Ri = 0 then halt else continue; 
set Ry := Ri. 


We can clock w- n by having n loops in separate lines of code, where loop! calls loop2 and loop 27 calls 
loop i+ 1 instead of the halting command, each time resetting Ri to 1. Instead, we could also use 
a separate register to perform a countdown from n. When trying to extend this approach trivially to 
clock w - w, we run into a problem that is connected to the fact that w -w is a compound limit, i.e. a 
limit of limits: At time w-w, Ri will have cofinally often been set from 0 to 1, so convergence or Ri 
fails. While this is easily fixed as seen below, it hints at the limitations imposed by the strict limit rule 
and why the supremum of ITBM clockable ordinals might be quite low in the ordinals. 


Lemma 2. The first compound limit ordinal w -w is ITBM clockable. 
Proof (by algorithm). 
Algorithm 3 


set Ri := 1; 

set Ro := 1; 

call loop; 

loop; 

if Ro = 0 then halt else continue; 

if Ri =0 then set Ro := ie and set Ri := Re and continue, else continue; 
set Ry := i; 

call loop; 


In this algorithm, Rı is halved repeatedly to detect limits. Once a limit time is reached (Ri = 0), Re 
is halved and Rj is reset not to 1 but to the current value of R2. Once Rə hits 0, we have found the 
compound limit w-w. At every limit, every register content converges. 


As before, it is easy to clock finite multiples of the form w-w-n. Also, if we extend the algorithm to 
use extra registers R3, R4,..., Rn in the same manner as we extended the w-algorithm with R2, we 
can in fact clock an finite power w”. However, this is as far as we can get, as w” turns out to be the 
supremum of the ITBM clockable ordinals. 


First observe that at any limit time, the register contents are a fixed point for every command that has 
been carried out cofinally often: 


Lemma 3. Let (Ci)i<o be an ITBM computation by some program P and let a < 6 be a limit ordinal. 
Then the first command in the computation that alters the register contents after time a has not been 
carried out cofinally often below a. 


Proof. Let (c,¢,0,:,-) € P be a computation node which is called cofinally often below a. Let c be 
called at some time 8 > a, where for all a < y < 8 the register content has not been changed yet, i.e. 
R(a) = R(y) = R(8). Since ¢ may be assumed as locally continuous (cf. Note 1), we get: 


R(B +1) = 4(R(3)) 
= 6(R(a)) 
= ¢(lim R(t)) 


tra 


lim (R(t) 


tr>aAI(t)=c 


= lim R(t+1) 


toaAI(t)=c 


= R(a) 


So the first computation node that changes the register content after time a cannot have been called 
cofinally below a. 


Theorem 1. Let P be a program with k computation nodes. Then in any computation (Ct)t<o according 
to P, the register contents stabilize before w*t?, 


Proof. If k = 0 then the computation halts after finitely many steps or diverges since the program 
contains only finitely many branch nodes: The computation may run through these nodes in finite 
sequence or in an infinite loop. Every branch node may trigger halting depending on its rational 
function evaluated on the unchanged input. After finitely many steps, every node in the sequence or 
loop has been visited once. If the computation didn’t halt up to this point, the program will go on 
forever as the register content is never changed. 


So let the hypothesis hold for k. Let P be a program with (k+1)-many computation nodes. Suppose the 
register contents change after w*t?. The first computation node c responsible for a new register content 
is not used cofinally often below w*t?. Let a be the supremum over the times < w*t? when c was carried 
out nontrivially. We can view (Ct)a<t<o as the a computation by P \ {c} on input R(a). Inductively, 
the register content of this computation stabilizes in w*t'-many steps. Since a + w*t! < wt? this 
means that the original computation stabilizes before w*t?. But a computation that stabilizes before 
a limit can never change its register content again. 
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Once the register contents have stabilized, an ITBM computation may run for an additional finitely 
many steps before halting, or diverges. So we get: 


Theorem 2. The supremum of ITBM clockable ordinals is w®. 
The above argument is in fact independent of the input: 
Corollary 1. Every ITBM computation halts before w“ -many steps or diverges. 


Corollary 2. If an ITBM computation diverges, the register content at time w” is not changed any 
more and can be considered the pseudo-output of the diverging computation. 


3 Connections to other models of computation 


From a computability perspective, reals provide ample possibilities as codes for complex objects. We 
can code and decode into reals with our ITBMs by interpreting the binary expansion of a real in the 
interval [0,1] as an element of the Cantor space “2, i.e. an w-long sequence of 0’s and 1’s. The binary 
expansion of a real is not necessarily unique, so two binary strings representing the same real will 
appear to our machines as equivalent. 


Let us give an algorithm to retrieve the n-th binary digit of an x € [0,1]. 


Algorithm 4 
input Ry = 2; 
input Ro =n; 
set Rg = 1; (= 2°) 
set R4 = 0; (the current approximation to x) 
call loop; 


loop: 

set Ro := Rə — 1; (= #remaining_ iterations) 

if Ra := 0 then call lastloop else continue; if Ra + Ra > x then continue (do not add Ra to the 
approximation R4) 

else set Ra := Ra + Rs and continue; (add Rz to the approximation) 

set Rs _ 4s (= DA iterations) 

call loop; 


lastloop: 
if Ra + Rs > z then set Ri := 0 and halt else set Ri := 1 and halt; 


With this algorithm we can also do local changes to the binary bits of x in the fashion of a Turing 
machine. So finite Turing computations can obviously be implemented on BSS/ITBMs. Also, the halting 
problem for Turing machines becomes ITBM computable: 


Lemma 4. The classical halting problem is ITBM computable. 


Proof. Since standard Turing computations are ITBM computable, we can generate Turing programs 
successively. So, in iteration n carry out the first n Turing programs for n many steps on empty input, 
using a dedicated simulation register. In step n, use only the n-th and later binary digits for the Turing 
simulation, so at time w, this register will have converged to 0 (cf. Note 2). Once the algorithm finds 
that some program (say, the i-th) halts on its input, change the i-th binary digit of an initially zero 
output register to 1. Since there is a finite time when all halting computations of programs with index 
< n will have halted, this register converges to the halting problem. 


So our machines have more computing strength than Turing machines or classical BSS machines. Also, 
the type-2 Turing machines of computable analysis (see [5]) can easily be simulated by ITBMs: The 
output tape of a type-2 Turing machine, when modelled as an ITBM register, converges by definition. 
Input tapes do not change their content and the finitely many work tapes can be made convergent like 
in Note 2. 


Using above coding, ITBMs can also operate on functions. A continuous function f : R —> R may be 
input to an ITBM as its restriction to the rationals f° = f | Q : Q > R coded into a real xy € [0,1]. 
This requires a fixed enumeration of rational numbers q : w — Q which may be chosen as ITBM 
computable and an ITBM computable bijective pairing function (-,-) : w x w —> w. Then f@(z) = y 
may be expressed as “the (q~+(a),q~+(y))-th binary digit of py is 1” (and 0 in the case of f(x) Æ y). 
By computing a convergent sequence of rationals for a given real (nested intervals) we can compute the 
function value at this real. This takes w-w-many steps: Produce, for all n < w, the approximations of x 
up to n binary digits in a similar way to Algorithm 4 as the desired sequence of rationals converging to 
x. For every such approximation, the function value needs to be decoded from py. This takes w-many 
steps each and can be done with bounded scratch registers. So, with the tricks described in Note 2, all 
registers can be arranged to converge at limits. 


Example 3. The derivative of a differentiable function is ITBM computable. 


Proof. Given a point x, have an ITBM evaluate the differential quotient in x using only f(x) itself and 
rational approximations to f(x). 


An upper bound on strength of ITBMs is given by the strength of ITTMs: 


Lemma 5. Let P be an n-register BSS program. There is an ITTM program Q and a map f : R” > “2 
s.t. for every x E€ R the ITTM computation by Q on input f(x) halts and returns the information that 
either no computation by P on x exists, or that P halts on x with output y € “2, or that the P diverges 
on x with pseudo-output y € “2, where f—'(y) is the final register content of the ITBM computation 
by P ona. 


Proof. We assume that the ITTM we are working with has a finite number of read-write tapes, which 
can be accomplished of interlacing these tapes onto the single scratch tape in the definition of ITTMs 
in [3]. Due to Corollary 2, we know that after time w” also a diverging ITBM computation does not 
change its register content anymore. Since it is easy for an ITTM to construct a well order of w of length 
w” on one tape, it is possible for an ITTM to code the complete ITBM computation (Rz, It)t<w»e up 
to time w” on the output tape. 


So let us begin with defining a map f’ : R > “2. First, imagine the value f'(x) to be contained in three 
elements of “2 (i.e. three Turing tapes) where the first tape contains only a 0 or 1 in the first cell to 
specify the sign of x, the second contains the maximal exponent n s.t. 2” < x, and the third contains 
the binary expansion of žr, ignoring the decimal (binary) point and normalized in the following way: 
Binary expansions that end on an infinite trail of 1s are replaced by the unique binary expansion of the 
same number that ends on a trail of 0s. Then, use an ITTM computable pairing function to interlace 


these three tapes into one. The function f’ can easily be extended to a function f : R” > “2. 


It is easy to see that an ITTM is perfectly capable - albeit with an enormous time consumption - 
of computing addition, subtraction, multiplication and division on elements of “2 and of normalizing 
the result in the way described above. So, given the input, the program Q will start carrying out the 
sequence of ITBM commands in P and writing the results (and the program instructions used) one 


after another in the respective cells of the output tape. At limit times, the output tape contains all the 
information of the previous times, so it is easy for the Q to check whether everything converges and 
compute the limit if one exists. If not, output that there is no valid computation by P on a. If yes, Q 
will continue its simulation of P up to time w”. At time w”, it can replace the output tape content 
with f(y) where y is P’s output or pseudo-output. 


3.1 Open questions 


1. 


Details in the definition of ITBMs might have strong influence on the strength of the resulting 
machines. Other (stronger) variants may have interesting properties as well. 


. When allowing parameters, it appears that functions higher up in the Baire hierarchy of iterated 


pointwise limits of continuous functions might be computable by our machines. 


. With previous models of ordinal computations, tie ins to descriptive set theory and constructibility 


theory proved to be very insightful. We believe this to be applicable to ITBMs as well. 
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